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1. Introduction 


In this study we present results of a series of numerical experiments 
in which an upper ocean model l.s driven by surface heat fluxes and stress 
fields derived from the FGGE SOP-1 GLAS analysis/forecast system (Halem et 
al. , 1982). 

The ocean model is essentially an oceanic boundary layer coupled with a 
dynamic model of the upper ocean currents (Schopf anu Cane, 1983) • The 
oceanic boundary layer can be viewed as an intermediary zone between the 
deep ocean circulation and the atmosphere, extending from the surface to a 
depth of about 10-150 meters. In this layer both the temperature and the 
salinity fields are almost constant in the vertical. This vertical constancy 
is due to mixing which is caused by turbulence. The turbulence is in turn 
caused by breaking waves, shear Instabilities, etc. 

The aim of this study is to seek a better understanding of the model 
ocean boundary layer response to various forms of atmospheric forcings and 
the associated time and space scales involved. It is known that the time 
scales of the atmospheric flow are smaller than those of the oceanic flow 
even for an upper ocean model like ours. On the other hand, the space scales 
for the baroclinic motion in the ocean are considerably less than in the 
atmosphere. The grid being used is of 4° by 5° of latitude and longitude, 
respectively and was chosen £o coincide with that of the atmospheric analysis. 
However, an accurate description of the oceanic synoptic scales requires a 
much finer horizontal resolution, on the order of 10-50 kilometers. The coarse- 
ness of the grid is very significant limitation of our experiments. 

Although the distributions of atmospheric wind forcings are very similar 
in summer and in winter, their intensity is greater in winter (Elsberry and 
Camp, 1978: Elsberry and Raney, 1978). The source of mechanical energy is 
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proportional to u^. . Thus, during the passage of a winter storm, there Is 
an upward heat flux and the mechanical mixing is enhanced by the strengthening 
of the " 7 ind forcing. As a result of the upward heat flux, the bottom of the 
mixed layer deepens. After the storm's passage, the cooling effect continues, 
due to entrainment. In summer, the winds are weaker. Consequently, the 
entrainment at the base of the oceanic boundary layer Is almost nonexistent. 

As a direct consequence of the stabilizing effect of the surface heating, the 
oceanic boundary layer becomes shallower. (Elsberry and Camp, 1978: Elsberry 
and Raney, 1978). 

One of the problems faced in this kind of study is the intialization of 
the ocean model. A solution may be available in the near future: according 

to Levoy (1981), it is expected that global SST distributions with higher 
resolution will be possible to get from satellite radiances. Also, accurate 
estimates of the surface wind field may be obtained using data from the 
radar Scatterometer (O'Brien, 1981). However, In the present study somewhat 
arbitrary initial conditions are used for the mixed layer depth and deep 
temperature fields. The initial currents were assumed to be zero and the 
initial SST's were obtained from the GLAS/temperature retrieval system for 
January 1979 (Susskind et al . , 1982). It is very possible that the results 
may have been dependent on such a choice. 

Another characteristic of these experiments is that we have used a one- 
way coupling mechanism, in which the atmospheric parameters derived from the 
GLAS analysis/forecast system were used to drive the ocean model. It should 
be noted that the atmospheric analysis was performed assuming climatological 
SST's, and that the coupling did not allow any feedback from the ocean model's 
predicted SST to the atmosphere. The lack of feedback may also be a serious 
deficiency of our experiments. 


- 2 - 


ORIGINAL PAGE \$ 

2. The ocean model and atmospheric analysis OF POOR QUALIT 

The mixed-layer ocean model used in this stndy is the one developed by 
Schopf and Cane (1983). Following Kasahara (1974), the model equations of 
momentum, temperature, continuity and the hydrostatic relation are as follows: 


3 (hV)/3t + V* (VhV) 4* 3(w e V)/3s - If 4- u/a tane] k x hV- 

N /V rV tsj 

- h { Vp 4- b Vz} 4- 3 t/3s + h F H (V), (1) 

fs> rsJ 

3(hT)/3t 4- V« (VhT) + 3(w e T)/3s - - 3Q/3s - 3D v /3s + h D H (T) (2) 

rsf 

3h/3t + V«(hV) + 3w e /3s « 0 (3) 

r*t 

3 P/3 s * bh. (4) 

where s is a generalized vertical coordinate. The bouyancy, b, is calculated 
via a linear expression: 

b * b(T) * gct(T - T r ). (5) 

Following Niiler and Krauss (1977) and Kira (1976) the equation of entrain- 
ment, w e , needed to clo' a the system is: 

w e H(w e ) t h 1 (b 1 - b e ) + q 2 - m e |AV| 2 } - 2 m g uj - 2 e^hj 

+ hiB 0 ll - m b H(B 0 )J + B^h). (6) 


For a more detailed explanation about the model and the symbolism being 
used, the reader is refered to the original paper (Schopf and Cane, 1983). 

In these experiments the model has been used in a global configuration 
with a resolution of 4° of latitude by 5° of longitude. While no slip con- 
ditions are used at the coastal boundaries, no flux of temperature is allowed 
through these same boundaries. The time step used is of three hours. 

Because salinity is not predicted in the model, and in order to avoid 
complications associated with ice generation, artificial boundaries have 
been arbitrarily placed at 70°N and 58°S. 
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3. Description of the experiments 

The parameters that drive the ocean model are: the surface heat budget, 
the wind stress and the friction velocity. The components of the surface 
heat budget include the latent and sensible heat fluxes and the incoming and 
outgoing radiative fluxes. Thus, the surface heat balance, Q, is computed 
as : 

Q • (S* - X*, - E v - S H )/ P (7) 

where Sy, X^, E v and represent the solar radiation at ground level, the 
long wave radiation at ground level, the latent heat, and the sensible heat 
flux, respectively. A reference sea density is represented by p. 

The surface v/ind stress, ^ , is computed as in the atmospheric model while 
the friction velocity, u x , is calculated as usual (Somraerville et al. , 1974). 

Both Q, u x and T from the 6 hour cycle global analysis produced by the GLAS 
analysis/forecast system (Halem et al . , 1987.). 

The geopotential height, horizontal wind, and relative humidity lire anal- 
yzed at mandatory pressure levels in the GLAS objective analysis scheme 
(Baker, 1983). The 6-hour model forecast supplies a first-guess for the 
above fields at 300 mb and sea level. The assimilation/forecast model is a 
fourth-order global atmospheric model which is based on an energy-conserving 
scheme with all horizontal differences computed with fourth-order accuracy 
(Kalnay-Rivas et al. , 1977; Kalnay-Rivas and Hoitsma, 1979). 

Five experiments to study the result of using the FGGE GLAS analysis/ 
forecast system to drive (without feedback) a global upper ocean model, were 
performed for the period 5 January to 9 February, 1979 (FGGE first special 
observing period). 

Preliminary experiments coupling the atmospheric forcing from the GLAS 
FGGE analysis showed unrealistic large changes in the sea surface temperature 
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(SST). In order to understand the causes of these unrealistic variations, a 
series of simpler idealized experiments were conducted. The main objectives 
in setting up these experiments were to study tint time evolution and seek a 
better understanding of the time scales involved in the development of >. ( sym- 
metrlc oceanic features. 

In the first three idealized experiments, the oceanic initial conditions 
were zonally averaged (Figures 3.1 and 3.2). In addition, in the first experi" 
went, all the external atmospheric forcings were obtained from the January 1979 
analysis and zonally averaged; i.e., the surface heat balance, the friction 
velocity, were assumed to be constant in time and independent of longitude 
(Figures 3.3, 3.4, and 3.5). 

In the second experiment, the wind stress and friction velocities were 
zonally averaged throughout the entire time integration. For the surface 
heat balance, however we used the January 1979 average value (Figure 3.6) as 
obtained from the GLAS analysis. 

In the third experiment, the only forcing parameter that was zonally 
averaged was the surface heat balance. On the other hand, for the other two 
atmospheric parameters — friction velocity and wind stress — their respective 
January 1979 average values were used during the five week integration. 

In the last two experiments (fourth and fifth) , the initial condition for 
the mixed layer depth is depicted in Figure 3.7. In the fourth experiment, 
all the atmospheric forcing parameters were the time average values corresponding 
to January 1979 (Figures 3.6, 3.9, 3.10, and 3.11). 

Finally, in the last experiment the instantaneous values of all the externally 
atmospheric forcings were used in the time integration. 

The initial SST conditions were the January time averaged SST obtained by 
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Susskind et al . (1982) and is shown in Fig. 3.8. For comparison with the 


numerical results we also present the averaged February SST (Fig. 3.12) and 
the difference between the two fields (Fig. 3.13). Also, for the sake of 
completeness, the climatological difference between February and January 
is depicted in Fig. (3.14). 

In all five experiments, the ocean model was started from rest. 


- 6 - 



3.1 


Zonally averaged ML depth 
Contour interval of 5 m 
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3.5 


Zonally averaged wind stress in the x direction 
Contour interval of 75 10“-*N/sn* 
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3.13 


Observed sea 
January 1979 


surface temperature difference between February and 
(Susskind et al. 1982) Contour interval of 0.5° C 


19 






cO 60 90 120 150 

iperature difference between February 
of 0.5° C. 


4* Results 


ORIGINAL PAGE 13 
OF POOR QUALITY 


The discussions of these experiments will be limited to the open ocean 
regimes. Thus, t will avoid the study of the especially complex circulations 
in coastal regions, were the lack of boundary currents in the initial conditions 
and the coarseness of the grid, moke the model particularly unrealistic. 

First Case : (zonally averaged initial conditions, zonally averaged surface 
heat flux, wind stress and friction velocity) 

The analysis of Figures 4.1.1 and 4.1,5 shows that most of the changes 
in the ML depth occurs in the first ten days. The feature is observed 

for the horizontal velocity (Figures 4,1,3, 4.1.4, 4,1.7 and 4.1.8). However, 
the changes in the temperature field do occur in a much longer time (Figures 
4.1.2 and 4.1.6). More specifically, no net change of temperature is observed 
after one week. 

Since the results show that the principal changes after one week in depth, 
aud velocity field are essentially retained after five weeks for the first three 
eases , to avoid redundancy we will limit our discussion to on examination and 
description of the fifth week's results. 

In Figure 4.1.6 we observe a decrease in temperature north of 10°S up to 
45°N. This can be attributed to two factors: the deficit of surface heat 

balance north of 10°N (Figure 3,3), combined with equatorial upwelling (notice 
the divergence implied by Sv/3y in Fig. 4.1.8). Conversely, the increase in 
temperature south of 10°S should be attributed to the surplus of surface heat 
balance; north of 50°N the winds and heat balance are very weak, and so is 
the change in temperature. 
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Weekly averaged value of the zonal velocity field, u, after one week 
for the first case. Contour interval of 75 mm/s 






4.1.4 Weekly averaged value of the meridional velocity field, v, after one 
week for the first case. Contour interval of 75 mm/s 
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4.1.6 Weekly averaged the ML temperature changes after five weeks for the 
first case . Contour interval of 1 ° C 
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The weak winds north of 30°S are responsible for the shallowing of the 
mixed layer at those latitudes especially at high latitudes in the northern 
hemisphere (Figures 3.5 and 4.1.5). The deepening of the alxed layer in the 
belt between 40°S-55°S is due to the westerly winds, which are stronger in 
the southern hemisphere. 

Second Case : (zonally averaged initial conditions, zonally averaged wind stress 
and friction velocity, time averaged, but non zonally averaged 
heat flux) 

As in the previous case most of the changes in the height and in the 
horizantal velocity of the oceanic boundary layer model take place in the 
very first days (Figures 4.2.1, 4.2.3, 4.2.4, 4.2.5, 4.2.7 and 4.2.8). Again, 
from Figure 4.2.2 it is observed that the temperature field changes observed 
after one week are small. 

There is a marked difference in the temperature field between this case 
and the previous one. From a comparison of Figures 3.6 and 4.2.6, we observe 
that whenever we have a deficit of surface heat balance, the ocean model 
responds by a decrease in the temperature field, and viceversea. The only 
exception is again north of 50° N. 

In the northern hemisphere, where we have a maximum deficit of surface 
heat balance, we observe a deepening of the mixed layer (Figure 4.2.5). There- 
fore the deepening of the mixed layer observed in the belt of 15°N-35°N in 
the sub-tropical Atlantic and Pacific Ocean is due to both the deficit of 
surface heat balance and the strength of the winds observed at those latitudes 
(Figures 3.5 and 3.6). 

The analysis of these same figures demonstrates that the effects of the 
lighter winds north of 35°N are more important than the effects of the deficit 
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Weekly averaged ML depth changes after five weeks for the second case 
Contour interval of 10 m 
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4.2.6 Weekly averaged ML temperature changes after five weeks for the second 
case. Contour interval of 1°C 
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of surface heat balance. This is reflected in the shallowing of the mixed 
layer observed north of 35°N (Figure 4.2.5). Similarly, the shallowing observed 
south of 15 °N should be attributed to the very moderate winds in this region 
and warming by the surface heat flux. As in the previous case, the shallowing 
of the mixed layer is interrupted by a deepening observed between 40°S and 55°S. 

It is also observed that the increase in temperature in that belt is somewhat 
Intense (2-3°C) and it can be stated that in the belt of 40°S to 55°S the effects 
of the surplus of surface heat balance and the strong winds do not tend to balance 
each other as far as the temperature field is concerned. 

Third Case : (zonally averaged initial conditions, zonally averaged heat flux, time 
averaged, but not zonally averaged wind stress and friction velocity) 

For the same reasons as before I will limit my analysis to the five week 
results. The results after one week (Figs. 4.3. 1-4. 3. 4) are also included 
for reference. 

In Figure 4.3.5, a deepening of the mixed layer In the north-central Pacific 
and extra-tropical northern Atlantic is observed. Thus the zonal asymmetry in 
the oceanic boundary layer height reflects the zonal asymmetry in the friction 
velocity (Figures 4.3.5 and 3.9). 

Since the deficit of solar radiation north of 40°N is weak (Figure 3.3), 
no net change in the deepening or shallowing of the mixed layer occurs. The 
shallowing of the mixed layer west of South and Central America is due to the 
presence of very light winds. The same effect is observed In most of the 
South Atlantic. In this particular case, in the belt between 40°S and 60°S, 
the deepening of the mixed layer is stronger than in the previous cases as a 
direct result of stronger westerly winds. The shallowing of the mixed layer 
obe .ved in the northern and central Indian Ocean is due to the presence of 
very light to moderate easterly winds, combined with a surplus of solar radiation. 
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4.3.2 


Weekly averaged ML temperature changes after one week for th© third 
case. Contour interval of 1° C 
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Weekly averaged value of the zonal velocity field, u, after one week 
for the third case. Contour interval of 75 mm/s 
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Weekly averaged ML depth changes after five weeks for the third 
case. Contour interval of 10 m 
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The sharp decrease in temperature in the Atlantic and Pacific Oceans 
at the equator should be attributed to the equatorial upwelling (Figures 4.3.6, 
4.3.7 and 4.3<,8) and by advection by the meridional wind. The increase in 
temperature south of 30°S should be to the surplus of the surface heat balance, 
the effects of which outweigh the effect of the wind stress. On the other hand, 
in the Atlantic and Pacific the decrease in temperature observed in the northern 
hemisphere up to 45°N should be attributed to the deficit of surface heat balance 
(Figure 3.3). 

Fourth Case : (Time averaged but not zonally averaged heat flux, wind stress 
and friction velocity) 

In these last two cases we will examine the differences between the first 
and fifth week results in more detail. 

The first set of pictures we are going to analyze details the mixed layer 
depth and temperature fields after one week, Figure 4.4.1 shows a uniform global 
deepening of the mixed layer depth in the latitudes 40°S to 50°S, At the same 
time and in the same place, we can observe a slight cooling (Figure 4.4.2). In 
this region the effects of the strong westerly winds are almost balanced by the 
surplus of surface heat balance. 

The shallowing of the mixed layer observed west of Central America and in 
the South Pacific is due to the combination of light winds and the surplus of 
surface heating. An increase in the SST is also observed. The deepening of the 
mixed layer observed in the central tropical Pacific can be explained by the 
stronger and more persistent easterlies (Figure 3.10). No net change of temper- 
ature is observed in that region, since the effects of the surplus of surface 
heat balance tend to counterbalance the effects of the stronger easterly winds. 

A shallowing of the mixed layer in the north and central Indian Ocean is 
caused by the effects of the light and variable winds, and the surplus of 
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Weekly averaged value of the meridional velocity field, 
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Weekly averaged ML depth changes after one week for the fourth 
case. Contour interval of 10 m 
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Weekly averaged ML temperature changes after one week for the 
fourth case. Contour interval of 1° C 






surface heat balance. These two factors also contribute to the slight increase 
in temperature observed in that region. Similar results can be found for 
most of the south Atlantic. 

An analysis of the second set of pictures indicates no demonstrable 
differences in the mixed layer depth between the first and fifth weeks (Figures 
4.4.1 and 4.4.5). The decrease in temperature observed in the equatorial Pacific 
is due to the strength of the local equatorial upwelling and the strong easterlies 
(Figures 4.4.6 and 4.4.8). Because of the surplus of surface heat balance, we 
observe a significant increase in the temperature field both in the south Atlantic 
and in the extra- tropical latitudes of the South Pacific, as well as along the 
western coast of South America. The results after one week are depicted in 
Figs. 4.4. 1-4. 4. 4; and after five weeks in 4. 4. 5-4. 4. 8. 

Because the effects of the westerlies are compensated by the surplus of 
solar heating, which is very strong at these latitudes, the decrease in temper- 
ature in the southern belt between 40°S and 50°S is marginal. The decrease in 
temperature observed both in the tropical Atlantic and Pacific Oceans in the 
northern hemisphere can be explained by the deficit of surface heat balance and 
the light and persistent easterlies. 

Fifth Case : (Instantaneous values for all atmospheric forcing) 

We will begin with an analysis of the one week results (Figs. 4.5. 1-4. 5.8) . 

In this case, in the belt between 35°S and 60°S the deepening of the mixed layer 
Is well organized globally, (Figures 4.5.5). The changes in temperature after 
one week are (Fig. 4.5.6) are generally smaller than 1°C. The main exception is 
near the East coast of North America, It is interesting to notice that the use 
of time varying wind stress and friction results in smaller changes in ML depth 
(Fig. 4.4.1). A discontinuous temperature front can be observed to the Southeast 
of Africa. This front is stronger than in the previous cases. 
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4.4.4 Weekly averaged value of the meridional velocity field, v, after 
one week for the fourth case. Contour interval of 75 mm/s 
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Weekly averaged value ML temperature changes after five weeks for 
the fourth case. Contour interval of 1° C 
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4.4.7 Weekly averaged value of the zonal velocity field, u, after five 
weeks for the fourth case. Contour interval of 75 mm/s 
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Weekly averaged ML temperature changes after one week for the fifth 
case. Contour interval of 1° C 




After five weeks the deepening at the ML observed in most of the south 
Atlantic may be explained by the deficit of surface heat balance which counter- 
balances the effects of the very light easterlies observed in the region. 
Opposite reasons are valid for the shallowing of the ML observed in the south 
Pacific down to 40°S (Figures 4.5,9 and 4.5.13). Also very light easterly 
winds are observed in the south Pacific after five weeks (Figures 4.5.11). 

As far as the temperature in the north and central Indian Ocean is concerned, 
the deficit of a strong surface heating balance, combined with the minimal 
balancing effects of the light winds, produces a decrease in the temperature 
fields. The light winds do not balance the deficit of a strong surface 
heating balance, thus decrease in the temperature fields (Figures 4.5.14). 
Because of the absence of easterly winds west of southern Africa, we observe 
a great increase in the SST. The same is true of the west of South America. 
Furthermore, north of 10°S the deficit of heat balance produces a decrease 
in the SST in the Atlantic Ocean up to 35°N. Results for the fifth week are 
depicted in Figs. 4.5.8-4.5.14. 

This decrease in temperature is not as remarkable in the Pacific Ocean 
because of a strong surplus of surface heat balance. The presence of moderate 
easterly winds tends to counterbalance the effects of the surplus of surface 
heat balance in the belt between 10°N and 30 6 N in the Central Pacific, explain- 
ing the slight decrease in the SST observed in that area. North of 30°N, in 
the same Pacific Ocean, the surplus of surface heating tends to counterbalance 
the effects of the westerlies, thus generating no net change in the temperature 
fields. 

Because the time scale of response of the mixed layer temperature to the 
atmospheric wind forcing is a matter of hours, the instantaneous effects of the 
stronger wind stress will have a greater impact than in the time averaged case. 
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4.5.10 Weekly averaged value of friction velocity for the fifth week 
Contour interval of 15 10“ 1 mm/s 
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4.5.11 


Weekly averaged value of the wind stress in the x direction for the 
fifth week. Contour interval of 75 lCT^N/m 























Then, after five weeks the temperature front observed in the Indian Ocean at 
35°S is stronger than after one week. Simultaneously, the decrease in tempera- 
ture observed after five weeks in the belt between 40°S and 60°S is even greater 
than after one week. East of Asia, the decrease in SST can be explained by 
the deficit of surface heating in that area. 

We now compare the five week temperature changes predicted by the model 
(Fig. 4.5.14) with the actual observed changes from January to February 1979, 
(Fig. 3.13). First of all we should point out the model has less skill in 
predicting the changes than climatology (Fig. 3.14). Although there is some 
skill in predicting the sign of the change, their magnitudes are generally 
overpredicted. The model has correctly predicted cooling in the North Atlantic 
and North Pacific. The observed warming west of South America and South Africa 
are also predicted, but larger than observed. The cooling in the South Indian 
Ocean and in the Australian region is well predicted, but the model predicted 
cooling in the Arabian Sea and Bay of Bengal, whereas warming was observed. 
Similarly, the model overpredicted the cooling northeast of South America. 

A major failure occured in the belt between 30°S to 60°S, where there was 
generally observed warming, with some small regions with cooling. Hare the 
model predicted strong cooling, presumably through excessive deepening of 
the mixed layer by the roaring forties. Results after five weeks are depicted 
in Figs. 4.5.15 and 4.5.16. 

It is interesting to note that, of all five experiments, the first two 
have the best agreement with the climatological changes (Fig. 4,1.6 and 4.2.6). 
The second one (in which surface friction and wind stress were zonally and 
time averaged and the surface heat flux was only time averaged) has the best 
skill in predicting the observed changes (cf Fig. 4.2.6 with Fig. 3.13), 
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Conclusions 


We have performed a series of experiments with the Schopf and Cane 
(1983) upper ocean model driven by surface fluxes derived from the GLAS atmos- 
pheric analysis for January 1979, with SST’s derived for January 1979 by 

Susskind et al . (1982), and starting from a state of rest. 

The model results are generally qualitatively reasonable. Namely, that 
whenever we have strong winds a deepening of the oceanic boundary layer is 
observed, and viceversa. Conversely, the opposite is also true. Furthermore, 
the temperature changes respond to the surface heat budget. Whenever an 
upward heat flux is observed, the boundary layer model responds by a decrease 
of temperature and viceversa. 

The model results show that most changes in the mixed layer height and 

horizontal velocity occur in the first days. On the other hand, changes in 

the temperature field take a longer time to develop. The best forecast was 
obtained in the experiment where the driving surface stress and friction 
velocity were zonally and time averaged, and the surface heat flux was time 
averaged. In the most realistic case (real initial conditions, instantaneous 
forcing fields from the atmospheric analysis) , the resulting changes in 
temperature were larger than observed and the correlation between observed 
and predicted changes was poor. 

The deficiency in the forecast of SST changes may be due to several 
factors: lack of sufficient ocean resolution, improper initialization, lack 
of feedback between the irsnsan and the atmosphere and the absence of transports 
by the strong boundary currents and perhaps unrealistic surface fluxes of 
heat and momentum. Unless these problems are alleviated it will not be 
reasonable to perform coupled atmospheric ocean forecasts. 
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Appendix: Computational implementation of the coupling 


The initialization of the model is in the sub-routine CLIMIN. The 

sub-routine reads ocean data from a data file. The fields are as follows: 

HI: Oceanic boundary layer depth; 

H2: Subsurface layer depth; 

Tl: Mixed layer temperature 

T2: Sublayer temperature 

T200: Temperature at the base of the subsurface layer; 

T300: Temperature at the pressure reference surface. 

The atmospheric forcing is fed into the model in the subroutine FORCIN. 
These values are taken from a log 8 tape. The values are as follows: 

TAU: Time in hours of data; 

TS: Ground temperature; 

U8, U9, V8, V9: Components of horizontal velocity at levels 8 and 9 in the 

GCM model; 

Ev: Evaporative heat flux; 

SH: Sensible heat flux; 

SW: Solar radiation fife ground level; 

XW: Long-wave radiation at ground level; 

WM: Wind magnitude; 

ST: Total surface wind stress. 

Also in this same subroutine the horizontal wind stresses are computed 
as TAUX and TAUY. The friction velocity and the surface heating budget are 
ftlso computed as USTR and HEATNG, respectively. 

The daily and weekly averages are computed in the main program. The 

variables are preceded by the letters PR. For example, the averages of 
i , H2, and Tl are FRU1, PRH2 f and PRT1 , respectively. 
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Figure Captions 


Figures 

3. 1 Zonally averaged ML depth 

Contour interval of 5 m 

3.2 Zonally averaged temperature 

Countour interval of 2°C 

3.3 Zonally averaged surface heat balance 

Contour interval of 20 Watts/in 2 

3.4 Zonally averaged friction velocity 

Contour interval of 15 lO^^mm/s 

3.5 Zonally averaged wind stress in the x direction 
Contour interval of 75 10“^N/m 2 

3.6 January 1979 average value of the surface heat balance 
Contour interval of 20 Watts/m 2 

3.7 Initial condition of ML depth (not zonally averaged) 

Contour Interval of 5 m 

3.8 Initial condition of ML temperature (not zonally averaged) 

(Susskind et al . 1982) Contour interval of 2° C 

3.9 January 1979 average value of the friction velocity 
Contour interval of 15 lO^^mm/s 

3.10 January 1979 average value of the wind stress in the x direction 
Contour interval of 75 lO^N/m 2 

3.11 January 1979 average value of the wind stress in the y direction 
Contour interval of 75 10"^N/m 2 

3.12 February 1979 average sea surface temperature (Susskind et al. 1982) 
Contour interval of 2° C 

3.13 Observed sea surface temperature difference between February and 
January 1979 (Susskind et al. 1982) Contour interval of 0.5° C 

3.14 Climatological sea surface temperature difference between February 
and January. Contour interval of 0.5° C. 

4.1.1 Weekly averaged ML depth changes after one week for the first case. 
Contour interval of 10 m 

4.1.2 Weekly averaged ML temperature changes after one week for the first 
case. Contour interval of 1° C 

4.1.3 Weekly averaged value of the zonal velocity field, u, after one week 
for the first case. Contour interval of 75 mm/s 
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4.1.4 Weekly averaged value of the meridional velocity field, v, after one 
week for the first case. Contour interval of 75 mra/s 

4.1.5 Weekly averaged ML depth changes after five weeks for the first case. 
Contour interval of 10 m 

4.1.6 Weekly averaged the ML temperature changes after five weeks for the 
first case . Contour interval of 1° C 

4.1.7 Weekly averaged value of the zonal velocity field, u, after five weeks 
for the first case. Contour interval of 75 mm/s 

4.1.8 Weekly averaged value of the meridional velocity field, v, after five weeks 
for the first case. Contour interval of 75 mm/s 

4.2.1 Weekly averaged ML depth changes after one week for the second case. 

Contour interval of 10 m 

4.2.2 Weekly averaged ML temperature changes after one week for the second 
case. Contour interval of 1°C 

4.2.3 Weekly averaged value of the zonal velocity field, u, after one week 
for the second case. Contour interval of 75 mm/s 

4.2.4 Weekly averaged value of the meridional velocity field, v, after one week 
for the second case. Contour interval of 75 mm/s 

4.2.5 Weekly averaged ML depth changes after five weeks for the second case. 
Contour interval of 10 m 

4.2.6 Weekly averaged ML temperature changes aftar five weeks for the second 
case. Contour interval of 1°C 

4.2.7 Weekly averaged value of the zonal velocity field, u, after five weeks 
for the second case. Contour interval of 75 mm/s 

4.2.8 Weekly averaged value of the meridional velocity field, v, after five weeks 
for the second case. Contour interval of 75 mm/s 

4.3.1 Weekly averaged ML depth changes after one week for the third case. 

Contour interval of 10 m 

4.3.2 Weekly averaged ML temperature changes after one week for the third 
case. Contour interval of 1° C 

4.3.3 Weekly averaged value of the zonal velocity field, u, after one week 
for the third case. Contour interval of 75 mm/s 

4.3.4 Weekly averaged value of the meridional velocity field, v, after one 
week for the third case. Contour interval of 75 mm/s 

4.3.5 Weekly averaged ML depth changes after five weeks for the third 
case. Contour interval of 10 m 

4.3.6 Weekly averaged ML temperature changes after five weeks for third 
case. Contour interval of 1° C 
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4.3.7 Weekly averaged value of the zonal velocity field, u, after five weeks 
for the third case. Contour interval of 75 mm/s 

4.3.8 Weekly averaged value of the meridional velocity field, v, after five weeks 
for the thivd case. Contour interval 75 mm/s 

4.4.1 Weekly averaged ML depth changes after one week for the fourth 
case. Contour interval of 10 m 

4.4.2 Weekly averaged ML temperature changes after one week for the 
fourth case. Contour interval of 1° C 

4.4.3 Weekly averaged value of the zonal velocity field, u, after one 
week for the fourth case. Contour Interval of 75 mm/s 

4.4.4 Weekly averaged value of the meridional velocity field, v, after 
one week for the fourth case. Contour interval of 75 mm/s 

4.4.5 Weekly averaged ML depth changes after five weeks for the fourth 
case. Contour interval of 10 m 

4.4.6 Weekly averaged value ML temperature changes after five weeks for 
the fourth case. Contour interval of 1° C 

4.4.7 Weekly averaged value of the zonal velocity field, u, after five 
weeks for the fourth case. Contour Interval of 75 mm/s 

4.4.8 Weekly averaged value for the meridional velocity field, v, after 
five weeks for the fourth case. Contour interval of 75 mm/s 

4.5.1 Weekly averaged value of the surface heat balance after one week. 

Contour interval of 20 Watts/m^ 

4.5.2 Weekly averaged value of friction velocity after one week. 

Contour interval of 15 10“ 1 mm/s 

4# 5.3 Weekly averaged value of the wind stress in the x direction after 
one week. Contour interval of 75 KT^N/m^ 

4.5.4 Weekly averaged value of the wind stress in the y direction after 
one week. Contour interval of 75 10 - 3N/m2 

4.5.5 Weekly averaged ML height changes after one week for the fifth 
case. Contour interval of 10 m 

4.5.6 Weekly averaged ML temperature changes after one week for the fifth 
case. Contour interval of 1° C 

4.5.7 Weekly averaged value of the zonal velocity, u, after one week for 
the fifth case. Contour interval of 75 mm/s 

4.5.8 Weekly averaged value of the meridional velocity, v, after one week 
for the fifth case. Contour interval of 75 mm/s 
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4.5.9 

4.5.10 

4.5.11 

4.5.12 

4.5.13 

4.5.14 

4.5.15 

4.5.16 


Weekly averaged value of surface heat balance for the fifth week 
Contour interval of 20 Watts/m^ 

Weekly averaged value of friction velocity for the fifth week 
Contour Interval of 15 10“^-nun/s 

Weekly averaged value of the wind stress in the x direction for the 
fifth week. Contour interval of 75 10"3N/m2 

Weekly average value of the wind stress in the y direction for the 
fifth week. Contour interval of 75 10”^N/m^ 

Weekly averaged ML height changes after five weeks for the fifth 
case. Contour interval of 10 ra 

Weekly averaged ML temperature changes after five weeks for the fifth 
case. Contour interval 1° C 

Weekly averaged value of the zonal velocity, u, after five weeks 
for the fifth case. Contour interval 75 mm/s 

Weekly averaged value of the meridional velocity, v, after five weeks 
for the fifth case. Contour interval 75 mm/s 
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